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Abstract 

We have considered a hot neutron star with a quark core, a mixed phase of quark-hadron matter, 
and a hadronic matter crust and have determined the equation of state of the hadronic phase and 
the quark phase, we have then found the equation of state of the mixed phase under the Gibbs 
conditions. Finally, we have computed the structure of hot neutron star with the quark core and 
compared our results with those of the neutron star without the quark core. For the quark matter 
calculations, we have used the MIT bag model in which the total energy of the system is considered 
as the kinetic energy of the particles plus a bag constant. For the hadronic matter calculations, 
we have used the lowest order constrained variational (LOCV) formalism. Our calculations show 
that the results for the maximum gravitational mass of the hot neutron star with the quark core 
are substantially different from those of without the quark core. 
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I. INTRODUCTION 



A hot neutron star is born following the gravitational collapse of the core of a massive star 
just after the supernova explosion. The interior temperature of a neutron star at its birth is of 
order 20 — 50 MeV Therefore, the high temperature of these stages cannot be neglected 
with respect to the Fermi temperature throughout the calculation of its structure. This 
shows that the equation of state of the hot dense matter is very important for investigating 
the structure of a newborn neutron star. Depending on the total number of nucleons, a 
newborn neutron star evolves either to a black hole or to a stable neutron star j^. Hence, 
calculation of the maximum mass of hot neutron star is of special interest in astrophysics. 

As we go from surface to the center of a neutron star, at sufficiently high densities, the 
matter is expected to undergo a transition from hadronic matter, where the quarks are 
confined inside the hadrons, to a state of deconfined quarks. Finally, there are up, down 
and strange quarks in the quark matter. Other quarks have high masses and do not appear 
in this state. Glendenning has shown that a proper construction of the hadron-quark phase 
transition inside the neutron stars implies the coexistence of nucleonic matter and quark 
matter over a finite range of the pressure. Therefore, a mixed hadron-quark phase exists 
in the neutron star and its energy is lower than those of the quark matter and nucleonic 

n 

matter p[. These show that we can consider a neutron star as composed of a hadronic matter 
layer, a mixed phase of quarks and hadrons and, in core, a quark matter. Recent Chandra 
observations also imply that the objects RX J185635-3754 and 3C58 could be neutron stars 
with the quark core j^. 

Burgio et al. have investigated the structure of neutron stars with the quark core at zero 
[sl and finite temperature js] with the Brueckner-Bethe-Goldstone formalism to determine 
the equation of state of the hadronic matter, they have used . We have calculated the 
structure properties of the cold neutron star by considering a quark phase at its core 3] 
and compared the results with our previous calculations for the neutron star without the 
quark core jsj. In these works, we have employed the lowest order constrained variational 
(LOCV) method for the hadronic matter calculations. In the present paper, we intend to 
extend these calculations for the hot neutron star with the quark core. 
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II. EQUATION OF STATE 

As it was mentioned in the previous section, we consider a neutron star composed of 
a hadronic matter (hadron phase), a mixed phase of quarks and hadrons, and a quark 
core (quark phase). Therefore, we should calculate the equation of state of these phases 
separately as follows. 



A. Hadron Phase 



For this phase of the neutron star matter, we consider the total energy per nucleon as 
the sum of contributions from the leptons and nucleons, 

E = Eiep + Enucl ■ (1) 

The contribution from the energy of leptons (electrons and muons) is 

Elep = Ee + E^, (2) 
where Ef. and are the energies of electrons and muons, respectively. 



E. = ^r vTT^ — 

n^nh^ Jo 1 + Exp{(3[miC^y/TTx^ - fii]} 

Here fii and rrii are the chemical potential and mass of particle i, /3 = (T is the 
temperature), n is the total number density of nucleons {n = Up + Un), c is speed of light 
and X is as follows, 

x= . 4 

niiC 

In our calculations, the equation of state of hot nucleonic matter is determined using the 
lowest order constrained variational (LOCV) method as follows |9l-ll6j . We adopt a trail 
wave function as 

^ = F^, (5) 

where (j) is the Slater determinant of the single-particle wave function and F is the correlation 
function which is taken to be 

F = 5n/(U)- (6) 

i>j 
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iS is a symmetrizing operator. For the energy of nucleonic matter, we consider up to the 
two-body term in the cluster expansion, 

Enud = Ei + E2. (7) 

The one body term Ei for the hot asymmetrical nucleonic matter that consists of Z 
protons and neutrons is simply the fermi gas kinetic energy, 

^1 = E ^« (8) 

i=l,2 

Label 1 and 2 are used instead of proton and neutron, respectively, and £i is 

E. = j:t^f^ik,T,n,), (9) 

^ ZTTli 



where f{k,T,ni) is the Fermi-Dirac distribution function 18| 



f{k,T,ni) ^p[,.(^k,T,ni)-,ii{T,n,)] ^ l' 

In the above equation, rii are the number densitis and are the single particle energies 
associated with the protons and neutrons, 

eAk.T.rii) = -, (11) 

where m* are the effective masses. 
The two-body energy, £'2, is 

^2 = ;^E<^^>(12)|U-J^>, (12) 



where 
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K12) = -|^[/(12), [V?2, /(12)]] + /(12)\/(12)/(12). (13) 



/(12) and V^(12) are the two-body correlation and inter-nucleonic potential. 

We note that the conditions of charge neutrality and beta stability impose the following 
constraints on the number densities and chemical potentials, 

rip = ne + rif, (14) 

l^n - f^p = fJ'e = f^fi- (15) 

The procedure to calculate the nucleonic matter has been fully discussed in the Refs. 



B. Quark Phase 



We use the MIT bag model for the quark matter calculations. In this model, the energy 
density is the kinetic energy of quarks plus a bag constant [B) which is interpreted as the 
difference between the energy densities of non interacting quarks and interacting ones jl?! . 

Stot = Su + Sd + Ss + B, (16) 
where £i is the kinetic energy per volume of particle z, 

= A r + h^k\y^f{k, T, n,)edk. (17) 

ZTT^ JO 

In above equation, g is the degeneracy number of the system and rij is the number density 
of particle i, 



n,; = A r f{k,T,n,)k^dk. (18) 
Jo 



27r2 JO 

For the quark phase, the Fermi-Dirac distribution function, f{k,T,ni), is given by 

^' ""'^ ^ Exp{(3i{mjc^ + h^^c^yl^ - /^.)} + 1 ■ ^^^^ 
We assume that the up and down quarks are massless, the strange quark has a mass equal to 
150 MeV and i3 = 90 MeVfm'^. Now, by applying the beta stability and charge neutrality 
conditions, we get the following relations for the chemical potentials and number densities, 

/Ud = /i« + /i«, (20) 

fls = f^u + IJ-i, (21) 

^ ^id = l^s, (22) 
2 11 

-Uu - -Ud - -Us -ni = 0, (23) 

ub = ^{nu + Ud + Us) , (24) 

where ni and /i; are the leptonic number density and chemical potential, and ri^ is the 
baryonic number density. 

The pressure of the system is calculated from free energy using the following equation, 

(25) 



where the Helmholtz free energy per volume (J-") is given by 

7 = Stot - TStof (26) 

The entropy of quark matter [Stot) can be written as follows, 

Stot = Su + Sd + Ss-, (27) 

where Si is the entropy of particle z, 

3 

Si{ni,T) = -kB [f{k,T,n,)ln{f{k,T,ni)) 

Jo 

+ (1 - f{k, T, n,)) ln(l - f{k, T, ni))]kHk. (28) 

C. Mixed phase 

For the mixed phase, where the fraction of space occupied by quark matter smoothly 
increases from zero to unity, we have a mixture of hadrons, quarks and electrons. In the 
mixed phase, according the Gibss equilibrium condition, the temperatures, pressures and 
chemical potentials of the hadron phase (H) and quark phase (Q) are equal j^. Here, for 
each temperature we let the pressure to be an independent variable. 

The Gibss conditions implies that 

= /^n , (29) 

^? = (30) 

where /i,^ and //^ (/i^ and yU^) are the neutron (proton) chemical potentials in the hadron 
phase and the quark phase, respectively, 

= d^J ^^'^ 

In above equations, S is the energy density of the system, 

£ = n{E + mc^). (33) 

To obtain and /j,^ for the hadronic matter in mixed phase, we use the semiempirical 
mass formula 
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E = T{n, x) + Vo{n) + (1 - 2x) ^2(71), (34) 



where a; = ^ is the proton fraction. T(n, x) is kinetic energy contribution and the functions 
Vq and V2 represent the interaction energy contributions which are determined from the 
energies of the symmetric nuclear matter {x = ^) and pure neutron matter {x = 0). We 
calculate Vq and V2 using our results for the LOCV calculation of nucleonic matter with the 
UVi4 + TNI nuclear potential which is discussed in section III A[ Now, we can obtain the 
chemical potentials of neutrons and protons from Eqs. (!3T|) - (1MI) as follows, 

/.f = T(n, x) + n^^^ + + Voin) + nV,'{n) 

+ (-3 + 8x - Ax'^)V2{n) + (1 - 2xfnV2\n) + mc^, (35) 



+ (1 - Ax^)V2{n) + (1 - 2xfnV2:{n) + mc^. (36) 
For the quark matter in mixed phase, we have 

/i^ = + /id, (37) 

/i^ = /i„ + 2iid- (38) 



At a certain pressure, we calculate yU„ for different /i^ under the condition that the densities 
yield this certain pressure. By calculating and /x^, we obtain /i^ and /i^. 

Now, we plot Hp versus /i„ for both hadron and quark phases, the cross point of the 
two curves satisfies the Gibss conditions. In the mixed phase, as the chemical potentials 
determine the densities, the volume fraction occupied by quark matter, x, can be obtained 
by the requirement of global charge neutrality, 

2 11 

Xi^nu - -nd - -Us) + (1 - x)np - rie = 0. (39) 
Finally, we can calculate the baryonic density of the mixed phase (M), 

nB = Xnq + (1 - (40) 
and then the total energy density of mixed phase is found, 

Sm = X^Q + (1 - X)^H. (41) 
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D. Results 



We have shown our results for the energy densities of hadron phase, quark phase and 
mixed phase in Figs. [T] and H] at two different temperatures. Figs. [T] and [2] show that at 
low densities the energy density of the hadronic matter is lower than those of other phases. 
However, as the density increases, at first the energy of mixed phase and finally the energy 
of quark phase is lower than those of other phases. We also see that there is a mixed phase 
for a range of densities. Below (beyond) this range, we have the pure hadron (quark) phase. 
By comparing Figs. [1] and [2l we see that for a given value of the density, the energies of all 
phases increases by increasing the temperature. 

Using the above calculated energy density, we can determine the equation of state and 
finally the structure of the hot neutron star with the quark core which is discussed in the 
next section. 



III. STRUCTURE OF THE HOT NEUTRON STAR WITH THE QUARK CORE 

The structure of neutron star is determined by numerically integrating the Tolman- 



Oppenheimer-Volkoff equation (TOV) 22|-|25| 



dP ^ G[£{r) + ^^][mir) + '-^] 

^ = A7rr^£ir), (43) 
ar 

where P is the pressure and S is the total energy density. For a given equation of state in 
the form P{S), the TOV equation yields the mass and radius of star as a function of the 
central mass density. 

In our calculations for the structure of hot neutron star with the quark core, we use the 
following equations of state: (i) Below the density of 0.05 fm^, we use the equation of state 



calculated by Baym 26|. (ii) From the density of 0.05 fm^ up to the density where the 
mixed phase starts, we use the equation of state of pure hadron phase calculated in section 
III A[ (iii) In the range of densities in which there is the mixed phase, we use the equation 
of state calculated in section HTCl (iv) Beyond the density of end point of the mixed phase, 
we use the equation of state of pure quark phase calculated in section III Bl All calculations 



8 



are done for B = 90 MeV fm~^ at two different temperatures T = 10 and 20 MeV. Our 
results are as follows. 

The gravitational mass as a function of the central mass density for the hot neutron star 
with the quark core at two different temperatures has been presented in Figs. E] and HI It 
is seen that for both relevant temperatures, the gravitational mass increases by increasing 
the central mass density and finally reaches a limiting value (maximum mass). In Figs. [3] 
and m our results for the case of neutron star without the quark core have been also given 
for comparison. We see that by including the quark core for the neutron star, our results 
for the gravitational mass are substantially affected. For the neutron star with the quark 
core, our results for the gravitational mass at three different temperatures (T = 0, 10 and 
20 MeV) have been compared in Fig. O It is seen that the gravitational mass increases by 
increasing the temperature. 

Figs. |6] and [7] show the gravitational mass versus the radius for both cases of the neutron 
star with and without the quark core at two different temperatures. At each temperature, 
it is seen that there is a reasonable difference between the mass-radius relations of these 
two cases of the neutron star. However, for both cases, we see that the radius decreases as 
the mass increases. By comparing Figs. [6] and U\ we can see that the decreasing rate of the 
radius versus the mass is substantially different for different temperatures. 

Our results for the maximum gravitational mass of the hot neutron star with the quark 
core and the corresponding values of radius and central mass density have been given in 
Tables [T] and [Tll for two different temperatures. Our results for the case of hot neutron star 
without quark core have been also presented for comparison. For different temperatures, it 
is seen that the inclusion of the quark core considerably reduces the maximum mass of the 
hot neutron star. This is due to the fact that by including the quark core for the neutron 
star, the equation of state becomes softer than that without the quark core. However, we do 
not see any substantial changes for the radius and central mass density of these two cases 
of the hot neutron star. 

IV. SUMMARY AND CONCLUSION 

For the hot neutron star, from the surface toward the center, we have considered a pure 
hadronic matter layer, a mixed phase of quarks and hadrons in a range of densities which 
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are determined by employing the Gibbs conditions, and a pure quark matter in the core, to 
calculate its equation of state at finite temperature. For calculating the equation of state of 
the hot hadronic matter, we have applied the lowest order constrained variational (LOCV) 
method at finite temperature. The equation of state of the hot quark matter has been 
computed using the MIT bag model with the bag constant B = 90 MeV fm^^ . Using this 
equation of state, we have solved the TOV equation by numerical method to determine the 
structure properties of the hot neutron star with the quark core at T = 10 and 20 MeV. 
Then, we have compared the results of these calculations with those for the hot neutron star 
without the quark core. It is found that our results for the maximum gravitational mass of 
the neutron star with a quark core are less than those of the neutron star without the quark 
core. 
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TABLE I: Maximum gravitational mass (Mmax)j and the corresponding radius (R) and central 
mass density (ec) of the hot neutron star without (NS) and with (NS+Q ) the quark core at 
T = 10 MeV. 

NS 2.07 10.22 26.94 
NS+Q 1.76 10.45 27.38 

TABLE II: As Table U but at T = 20 MeV. 
M^,,(Mo) R(Km) e,{W^^gr/cm^) 

NS 2.09 10.64 27.01 

NS+Q 1.78 11 27.37 
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FIG. 1: Energy density versus the baryonic density at T = 10 MeV 
for the hadron phase (solid curve), quark phase (dotted curve) and mixed phase (dashed curve). 
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FIG. 3: Gravitational mass versus the central mass density for the neutron star with (dotted curve) 
and neutron star without (solid curve) the quark core at T = 10 MeV. 
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FIG. 4: As Fig. Hbut at T = 20 MeV. 
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FIG. 5: Gravitational mass versus the central mass density for the neutron star with the quark 
core at r = (dotted-dashed curve), 10 (solid curve) and 20 MeV (dotted curve). 
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FIG. 6: Mass-radius relation for the neutron star with (dotted curve) and without (sohd curve) 
the quark core at T = 10 MeV. 
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